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ABSTRACT 

Aims. We present results from deep observations of the Galactic shell-type supernova remnant (SNR) RX J1713. 7-3946 (also known as 



O G347.3-0.5) conducted with the complete H.E.S.S. array in 2004. 



Methods. Detailed morphological and spatially resolved spectral studies reveal the very-high-energy (VHE - Energies £ > 100 GeV) gamma- 



^2 ray aspects of this object with unprecedented precision. Since this is the first in-depth analysis of an extended VHE gamma-ray source, we 

present a thorough discussion of our methodology and investigations of possible sources of systematic errors. 
^ Results. Gamma rays are detected throughout the whole SNR. The emission is found to resemble a shell structure with increased fluxes from the 
western and northwestern parts. The differential gamma-ray spectrum of the whole SNR is measured over more than two orders of magnitude, 
r\ , from 190 GeV to 40 TeV, and is rather hard with indications for a deviation from a pure power law at high energies. Spectra have also been 
■ determined for spatially separated regions of RX J1713. 7-3946. The flux values vary by more than a factor of two, but no significant change in 
■ " " ' spectral shape is found. There is a striking correlation between the X-ray and the gamma-ray image. Radial profiles in both wavelength regimes 
reveal the same shape almost everywhere in the region of the SNR. 

Conclusions. The VHE gamma-ray emission of RX J1713. 7-3946 is phenomenologically discussed for two scenarios, one where the gamma 
rays are produced by VHE electrons via Inverse Compton scattering and one where the gamma rays are due to neutral pion decay from proton- 
proton interactions. In conjunction with multi-wavelength considerations, the latter case is favoured. However, no decisive conclusions can yet 
be drawn regarding the parent particle population dominantly responsible for the gamma-ray emission from RX JI7I3. 7-3946. 

Key words, acceleration of particles - cosmic rays ~ gamma rays: observations - supernova remnants - gamma rays: individual objects: 
RX J1713.7-3946 (G347.3-0.5) 



1. Introduction 

It is commonly believed that the only sources capable of sup- 
plying enough energy output to power the flux of Galactic 
cosmic rays are supernova explosions (e.g., Ginzburg & 
Syrovatskii 1964). At the present time there are two main ar- 
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guments for this hypothesis: firstly, estimates of the power re- 
quired to sustain the observed nuclear Galactic cosmic -ray pop- 
ulation show that about 10% of the mechanical energy released 
by the population of Galactic supernovae would suffice, or, 
in other words, that supernova remnants could be the sources 
of the Galactic cosmic rays if the average acceleration effi- 
ciency in a remnant is about 10%. Secondly, a rather well de- 
veloped theoretical framework for the acceleration mechanism. 
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diffusive shock acceleration (for reviews see eg, Blandford & 
Eichler 1987; Jones & ElUson 1991; Malkov & O'C Drury 
2001), exists and it indeed predicts acceleration efficiencies in 
excess of 10%. 

The best way of proving unequivocally the existence of 
very-high-energy (VHE) particles, electrons or hadrons, in the 
shells of supernova remnants (SNRs) is the detection of VHE 
(about 100 GeV up to a few tens of TeV) gamma rays pro- 
duced either via Inverse Compton (IC) scattering of VHE elec- 
trons off ambient photons or in interactions of nucleonic cosmic 
rays with ambient matter. As was argued already in Drury et al. 
(1994), a system of imaging atmospheric Cherenkov telescopes 
with a large field of view provides the most powerful measure- 
ment technique for extended nearby SNRs at these very high 
energies. One should note that there exist two other experimen- 
tal approaches to trace VHE cosmic rays, the detection of X- 
rays, which suggests the presence of VHE electrons (Koyama 
et al. 1995), and of high-energy neutrinos, which probe exclu- 
sively nuclear particles. 

A prime candidate for gamma-ray observations is the SNR 
RX J1713. 7-3946, in particular because of its close associa- 
tion with dense molecular clouds along the line of sight (Fukui 
et al. 2003; Moriguchi et al. 2005), which might suggest 
a scenario of a supernova shell overtaking dense molecular 
clouds, leading to a detectable VHE gamma-ray signal from 
hadronic interactions, as described in Aharonian et al. (1994). 
RX J1713. 7-3946, situated in the Galactic plane, constellation 
Scorpius, was discovered in soft X-rays in 1996 in the ROSAT 
all-sky survey (Pfeffermaim & Aschenbach 1996). It is roughly 
70' in diameter and exhibits bright X-ray emission dominantly 
from its western shell. ASCA observations revealed that the 
X-ray emission is a pure non-thermal continuum without de- 
tectable line emission (Koyama et al. 1997; Slane et al. 1999). 
X-ray observations have also been conducted with Chandra 
and XMM with their superior angular resolution. Chandra ob- 
served a small region in the bright northwestern part of the 
SNR (Uchiyama et al. 2003; Lazendic et al. 2004). Despite 
distinct brightness variations within this small field, the corre- 
sponding X-ray spectra were all found to be well described by 
power-law models with similar absorbing column densities and 
pholon indices, albeit with rather large statistical uncertainties. 
XMM covered the remnant almost completely in five point- 
ings (Cassam-Chenai et al. 2004; Hiraga et al. 2005). Also on 
this much larger scale, a highly inhomogeneous and complex 
morphology was found in the western part of the SNR with two 
narrow rims resembling a double-shell structure running from 
north to south. The (non-thermal) X-ray spectra, when fit with 
a power law, exhibit strong variations in photon index across 
the remnant (1.8 < E < 2.6) and the hydrogen column den- 
sity A'h was found to vary significantly (0.4 X 10^^ cm~^ < 
A^H < 1.1 X 10^^ cm"^). The spectra of the central and the 
western parts differ clearly at low energies, possibly indicat- 
ing an increase in colunm density of AA^h ~ 0.4 x 10^^ cm"^ 
towards the west. Furthermore, a positive correlation between 
X-ray brightness and absorption was interpreted as being due 
to the shock front of RX J1713.7-3946 impacting a molecu- 
lar cloud in the west which was assumed to be responsible for 
the absorption. Further support for this scenario is lent by CO 



line emission observations with the NANTEN telescope (Fukui 
et al. 2003; Moriguchi et al. 2005), which suggest that the SNR 
is interacting with molecular clouds in this region at a distance 
of 1 kpc from the Solar System. The non-thermal X-ray emis- 
sion is possibly associated with interactions between the cloud 
and the western part of the SNR shell. 

Age and distance of the SNR are under debate and have 
been revised quite a few times. Initially, Koyama et al. (1997) 
had derived a distance of 1 kpc and correspondingly an age 
of about 1000 years from the colunm density towards the 
source as estimated from ASCA X-ray observations. Slane 
et al. (1999) on the other hand have derived a larger dis- 
tance of 6 kpc (corresponding to an age of about 10 000 years) 
based on the possible association of RX J1713.7-3946 with a 
molecular cloud in this region and the Hn region G347.6H-0.2 
to its northwest. Both the latest XMM and NANTEN find- 
ings are consistent with the remnant being closer, at 1 kpc, 
which might support the hypothesis of Wang et al. (1997), that 
RX J1713.7-3946is the remnant of a AD393 guest star which, 
according to historical records, appeared in the tail of constella- 
tion Scorpius, close to the actual position of RX J17 13.7-3946. 
The high surface brightness of this object, both in VHE gamma 
rays and non-thermal X-rays, suggests that it is close to the 
evolutionary phase where the shocks are most powerful. While 
hardly a conclusive argument, this implies that the remnant is 
observed at the sweep-up time when the ejecta are interacting 
with approximately their own mass of swept-up ambient mate- 
rial and the energy flux through the shocks (both forward and 
reverse) peaks. NormaUy this would be at an age of a few hun- 
dred to a thousand years, which indeed supports the closer dis- 
tance estimate. 

The radio emission of RX J1713.7-3946 is very 
faint (Lazendic et al. 2004) which puts it into a pecuUar class 
of shell-type SNRs with dominantly non-thermal X-ray and 
only very faint radio emission. The only other known object 
of this type is RX J0852.0-4622 (G266.2-1.2) (Aschenbach 
1998; Slane et al. 2001). 

RX J 17 13.7-3946 was detected in VHE gamma rays by the 
CANGAROO collaboration in 1998 (Muraishi et al. 2000) and 
re-observed by CANGAROO-II in 2000 and 2001 (Enomoto 
et al. 2002). Recently H.E.S.S., a new array of imaging at- 
mospheric Cherenkov telescopes operating in Namibia, has 
confirmed the detection (Aharonian et al. 2004b). This was 
the first independent confirmation of VHE gamma-ray emis- 
sion from an SNR shell. Furthermore, the H.E.S.S. measure- 
ment provided the first ever resolved gamma-ray image at very 
high energies. The complex morphology of RX J1713.7-3946 
was clearly unraveled. Together with the H.E.S.S. detection of 
RX J0852.0-4622 (Aharonian et al. 2005a) there are currently 
two spatially resolved VHE gamma-ray SNRs with a shell-like 
structure which agrees well with that seen in X-rays. These two 
objects may well be the brightest SNRs in the VHE gamma- 
ray domain in the whole sky; anything equally bright in the 
Northern sky would have been clearly seen in the Milagro sur- 
vey (Atkins et al. 2004), and the H.E.S.S. Galactic plane sur- 
vey (Aharonian et al. 2005d) reveals no SNRs brighter than 
RX J1713.7-3946 or RX J0852.0-4622 in the region covered. 
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The interpretation of the gamma-ray emission mechanisms 
for RX J1713. 7-3946 has been the subject of debate. From 
their flux level, the CANGAROO collaboration interpreted 
in Muraishi et al. (2000) the gamma rays as IC emission, 
whereas in Enomoto et al. (2002), after re-observations with 
CANGAROO-II, neutral pion decay was put forward as an 
explanation instead. The proposed model was then heavily 
disputed by Reimer & Pohl (2002) and Butt et al. (2002) 
because of its conflict in the GeV regime with the flux of 
flie nearby EGRET source 3EG 1714-3857 (Hartman et al. 
1999). Further attempts to model the broadband spectrum of 
RX J1713.7-3946 were undertaken (e.g., Ellison et al. 2001; 
Uchiyama et al. 2003; Pannuti et al. 2003; Lazendic et al. 
2004). However, they did not result in unequivocal conclusions 
concerning the acceleration mechanisms of the highest-energy 
particles or the origin of the VHE gamma rays from this source. 

Here we report on follow-up observations of 
RX J1713.7-3946 with the complete H.E.S.S. telescope 
array, conducted in 2004. The large field of view together with 
the high sensitivity of the system enable us to undertake for 
the first time detailed morphological studies in VHE gamma 
rays and measure spectral parameters in different regions of 
flie SNR. 

The paper is organised as follows: In Sect. 2 we present the 
data set and illustrate the performance of H.E.S.S. for obser- 
vations of extended gamma-ray sources. We explain in detail 
the analysis methods applied here in order to extract images 
and spectra. In Sect. 3 the results of spectral and morphologi- 
cal studies are presented along with systematic tests that have 
been performed in order to assure the validity of the analysis. 
Section 4 presents multi-wavelength data of RX J1713. 7-3946 
and its suiToundings aiming at putting our measurement at TeV 
energies into the context of the available data as preparation 
for a broadband modelling of the spectral energy distribution 
(SED) in Sect. 5. We discuss two scenarios for the generation 
of VHE gamma rays, a purely electronic and a purely hadronic 
one. The results are finally summarised in Sect. 6. 

2. Data processing 

2.1. H.E.S.S. observations 

Observations of RX J1713. 7-3946 were conducted between 
April and July 2004 with the High Energy Stereoscopic 
System (H.E.S.S.), a system of four imaging atmospheric 
Cherenkov telescopes (Hinton 2004; Hofmann 2003) situated 
in the Khomas Highland of Namibia, at 23° 16' S 16°30' E, 
1800 m above sea level. Each of the 13-m-diameter tele- 
scopes (Bernlohr et al. 2003; Cornils et al. 2003) has a tes- 
sellated Davies-Cotton mirror of 107 m^ area and is equipped 
with a 960-photomultiplier-tube camera (Vincent et al. 2003) 
covering a large field of view of 5° diameter. During stereo- 
scopic observations, an array-level hardware trigger requires 
each shower to be detected by at least two telescopes simulta- 
neously allowing for efficient suppression of the vast number 
of hadronic and muonic background events (Funk et al. 2004). 
The point source sensitivity reaches 1 % of the flux of the Crab 
nebula for long exposures («j 25 hours). 
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Fig. 1. X-ray image of RX J1713. 7-3946 (colour scale and thin grey 
contour lines, 1-3 keV, from Uchiyama et al. (2002)). The superim- 
posed thick white contours indicate the 94% and 98% levels of the 
detection-efficiency weighted H.E.S.S. exposure, given by the product 
of relative detection efficiency and the observation time. It can be seen 
that the relative gamma-ray detection efficiency between the centre 
region and the edges of the SNR differs only by x 5%. 



The observations were mostly performed in wobble mode 
around the SNR centre. In this mode the telescopes were po- 
sitioned such that the centre of the SNR was offset +0.7° in 
declination or right ascension away from the pointing direction 
of the telescope system, changing to the next position every 28 
minutes. Towards the end of the observation campaign, pure 
on-source pointings in which the centre of the SNR was co- 
incident with the system centre were additionally performed. 
In each of the five pointings the SNR RX J1713.7-3946, 
roughly 1° in diameter, was fully contained in the w 5° field 
of view of the system. The resulting effective exposure distri- 
bution (the product of the detection efficiency and the expo- 
sure) is illustrated in Fig. 1, where the ASCA X-ray measure- 
ment (Uchiyama et al. 2002) is shown (colour scale and thin 
contour lines) with superimposed white, almost circular con- 
tours indicating the 94% and 98% levels of the effective expo- 
sure. The observation strategy for this data set combined with 
the detector efficiency results in a very flat plateau in the re- 
gion of the SNR; from the centre to the boundaries the relative 
gamma-ray detection efficiency decreases by only about 5%, 
which is a great advantage compared to ASCA, for example. 
Not only is the SNR fully contained in all of the five point- 
ing positions, but one can also disregard for most purposes the 
modest change in relative detection efficiency from one region 
of the SNR to another. 
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2.2. Data sample 

The data comprise a total exposure time of 40 hours. Rejecting 
data taken under bad weather conditions, 36 hours of observa- 
tion time corresponding to 33 hours of Hve time remain for 
the analysis. The zenith angle of observations ranged from 
16° to 56° with a mean of 26°; it should be noted that about 
68% of the data were taken at small zenith angles between 
16° to 26°. The energy threshold (defined by the peak gamma- 
ray detection rate for a given source spectrum after all gamma- 
ray selection cuts) of the system increases with zenith angle. 
For the observations presented here, assuming a spectrum ap- 
propriate for RX J1713.7-3946, the threshold was « 180 GeV 
at 16°, «i 340 GeV at 40°, and 840 GeV at 56°. 

2.3. Data preprocessing 

The analysis technique applied here is described in de- 
tail in Aharonian et al. (2005b). After calibration of the 
data (Aharonian et al. 2004a), tail-cuts image cleaning is ap- 
pUed and the shower images in each telescope are parametrised 
in terms of their centre of gravity and second moments (HiUas 
1985). Stereoscopic event reconstruction based on the intersec- 
tion of image axes yields the shower direction, providing a res- 
olution of « 0.1° for individual gamma rays. Given the geome- 
try of the shower, cuts (optimised on Monte Carlo gamma-ray 
simulations and OFF source data, i.e. data without gamma-ray 
signal) are applied to select gamma-ray candidates and to sup- 
press the vast hadronic background. The gamma-ray energy 
is estimated from the image intensity and the reconstructed 
shower geometry yielding a resolution of « 15%. There is one 
main difference between a point-source and an extended-source 
analysis. In the latter case, the cut on the squared distance of 
events to the assumed source location is greatly increased in 
order to reflect the large source extension. In that case the en- 
ergy resolution worsens slightly (to $ 20% for a source of the 
size of RX J1713.7-3946) since an increased number of badly 
reconstructed events are included in the analysis. 

2.4. H.E.S.S. performance for extended sources 

The large field of view of the H.E.S.S. telescope system of x 5° 
diameter provides reasonable sensitivity for point sources at an 
angular distance up to 2° from the pointing direction of the tele- 
scope system (the point-source off-axis sensitivity derived from 
Monte Carlo simulations has been confirmed via observations 
of the Crab nebula, see Aharonian et al. (2005c)). Given the 
source diameter of up to 1 .3° and the offsets of 0.7° between the 
centre of the SNR and the telescope pointing direction during 
observations, it is important that the gamma-ray point-spread 
function is well behaved and does not broaden significantly 
with increasing offset from the pointing direction. Figure 2 
(left) shows the squared angular difference ^ between the re- 
constructed and the true direction of a simulated point source. 
The initial simulations have been generated at a number of 
fixed zenith angles between 0° and 70°. Taking then the zenith 
angle distribution of the actual data set and forming from that 
the weighted sum of the simulated distributions one obtains 



a resolution (taken as the 68% containment radius) of 0.08°. 
This is an order of magnitude smaller than the source diame- 
ter, which implies that in terms of angular resolution H.E.S.S. 
is well suited for morphological studies of RX J1713.7-3946. 
The right hand side of Fig. 2 illustrates the dependence of the 
angular resolution on the offset between source position and 
pointing direction for two representative zenith angles of 20° 
and 50°, for two different cuts on the minimum camera image 
size of 80 and 200 photo-electrons. It can be seen that the appli- 
cation of a higher cut on the image size improves the direction 
reconstruction by about 20% since only well defined camera 
images are used in the shower reconstruction, reducing fluctu- 
ation effects. However, the improved resolution is achieved at 
the expense of an increased energy threshold. In any case, for 
offsets smaller than 1.35°, which is the maximum offset under 
which parts of the SNR were observed, the resolution changes 
only slowly with offset and worsens by less than 10%. It was 
furthermore demonstrated with observations of the point-like 
source PKS 2155-304 (Aharonian et al. 2005b), that simulated 
point-spread functions agree well with data distributions and 
are thus well understood. 

The effective gamma-ray detection area depends on trigger 
conditions and analysis cuts. Well above the trigger threshold 
of the system, it is of the order of the area of the Cherenkov 
light pool (emitted by the secondary particle shower) on the 
ground. Typical effective area curves as a function of the offset 
angle between the gamma-ray source and the pointing direction 
of the system are shown in Fig. 3 for a zenith angle of 20°. The 
effective areas were determined from Monte Carlo simulations 
of a gamma-ray point source, accounting for the large extension 
of RX J 17 13. 7- 3 946 by increasing the cut on the maximum 
distance to the assumed source position. Since for an extended 
source like RX J1713. 7-3946 the flux has to be integrated over 
a larger solid angle than for point sources, the sensitivity of the 
system is reduced due to an increased background level. For 
moderate zenith angles and an integration region of 0.65° ra- 
dius around a gamma-ray source at 1° offset from the pointing 
direction (matched to the data set described here), simulations 
reveal that, for a source like RX J17 13.7-3946, the sensitivity 
drops roughly by a factor of four as compared to a point-source 
analysis. 

2.5. Analysis details 

The following section explains briefly the analysis applied to 
extract gamma-ray images and spectra of RX J1713.7-3946 
from the data. 

2.5.1. Morphology 

The gamma-ray images shown throughout the paper represent, 
unless otherwise stated, gamma-ray excess counts with back- 
ground subtracted. For the generation of these images a cut 
of 200 photo-electrons on the minimum camera image size in 
each telescope is applied to select a subset of events with su- 
perior angular resolution (see also Fig. 2 and Sect. 2.4). The 
background was estimated using about 160 hours of H.E.S.S. 
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Fig. 2. Left hand side: Squared angular distance between reconstructed and true direction for a Monte Carlo gamma-ray point source (point- 
spread function). A cut on the minimum size of images of 200 photo-electrons was applied. The simulated gamma rays followed a power 
law in energy with a photon index of 2. The zenith angle distribution and offset of the source with respect to the telescope optical axis was 
matched to the actual data set of RX J1713. 7-3946. The solid circles are the Monte Carlo histogram, the fit of a double Gaussian (dashed line) 
describes the point-spread function reasonably well. Indicated is the 68% containment radius of 0.08° (obtained from the fit), which is taken as 
the resolution of the data set. Right hand side: 68% containment radii of the point-spread function determined from simulations as function of 
the offset between the source direction and telescope pointing. Shown are simulations for 50° (triangles) and 20° zenith angle (circles) for two 
different size cuts (80 and 200 photo-electrons). The maximum offset of the actual data set of 1.35°, given by the wobble displacement of 0.7° 
and the approximate radius of the SNR of 0.65°, is indicated by the arrow. 
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Fig. 3. Effective detection area A for a moderate zenith angle of 20° 
as a function of offset angle 9 between the source location and the 
pointing direction of the system. Plotted are values at three represen- 
tative energy values of 0.2 TeV (near the post-cut trigger threshold), 
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simulated offsets effective area values are interpolated linearly. 



observations without any or only very faint gamma-ray sources 
in the field of view. All the events in these reference obser- 
vations passing gamma-ray cuts are assumed to be gamma-ray 
like background events and are used to estimate the background 
for the given data set. For that purpose the set of OFF runs has 
been divided into distinct zenith angle bands to account for the 
dependance of the system's gamma-ray acceptance on observa- 
tion altitude. Besides that, the acceptance depends only on the 
angular distance between shower and pointing direction and is 
to a very good approximation radially symmetric with respect 
to the pointing direction. Therefore, a radial 1-D lookup (num- 
ber of background events as a function of squared distance to 
the pointing centre) can be used in each zenith angle band for 
the background estimation. Given an observation at a certain 
zenith angle, a model background is created by selecting the 
1-D radial lookup from the zenith-angle band that matches the 
zenith angle of the observation. A 2-D background map of the 
sky is then created by rotating the coiTesponding 1-D lookup. 
Finally, the overall background map is created as the exposure 
weighted sum of the individual maps. A global normalisation 
factor a is calculated as ratio of the number of events in the data 
to the number of events in the background model, excluding 
regions that emit gamma rays. Figure 4 illustrates the validity 
of the approach for two H.E.S.S. data sets, RX J1713.7-3946 
(upper panel) and PKS 2155-304 (lower panel), the latter be- 
ing a point source for H.E.S.S.. Shown are slices along right 
ascension through the source centres. Overlaid on both data 
curves are the normalised background models of the whole 
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Fig. 4. As illustration of the good match between data and back- 
ground model, the data set of RX J1713. 7-3946 (upper panel) is com- 
pared to the high statistics data set of PKS 2155-304 (lower panel) 
which appears as a point source for H.E.S.S.. Shown are 1.3° and 0.2° 
slices along the RA axis through the centre of RX J1713. 7-3946 and 
PKS 2155-304, respectively, for the data (black) and the background 
model (green). The good match both in shape and absolute level be- 
tween the two curves illustrates the validity of the approach. The steps 
in the distributions (in the case of RX J1713.7-3946 at ±1.1° and 
+2.1°) are artefacts of the analysis: the usable range in the field of 
view of every observation was restricted to a radius of 2° (out of 2.5° 
physical camera radius) around the camera centre to allow only for 
well reconstructed events with good camera acceptance, and the fig- 
ure combines data from different pointings. 



data sets. In both cases, at different regions in the sky, for an 
extended and a point-like gamma-ray source, there is clearly 
a good match between model and data in regions outside the 
gamma-ray sources. 

Images of the gamma-ray excess are obtained by subtract- 
ing the normalised background model from the data in each bin 




-42 



17h20m 17h10m 17h00m 



Fig. 5. Count map of gamma-ray candidates for the region around 
RX J1713. 7-3946. A size cut of 200 photoelectrons on the camera 
images was applied. The bins are uncorrelated and the background is 
not subtracted. The white dashed circle indicates the region used to ex- 
tract the spectrum of the whole SNR (the ON region), the red circles 
indicate the OFF regions, the regions used for background estimation 
for the spectrum in each of the four wobble observation positions, 
which are marked as white crosses. The dashed red circle to the north 
of the SNR indicates an OFF region that was not used in the spec- 
tral analysis because it contains a gamma-ray source discovered in the 
H.E.S.S. Galactic plane survey (Aharonian et al. 2005d). 

/ of the 2-D map: 

^excess,! — ^data.i ^ ^background,! ■ 

Subsequently, these images are smoothed with a Gaussian to 
reduce statistical fluctuations. 

Typically the standard deviation for the smoothing is 
matched to the resolution of the data set, namely 30% - 50% 
of the 68% containment radius of the point-spread function. 
Here, for the data of RX J1713. 7-3946, a smoothing radius of 
2' is used. The resulting count maps are in units of integrated 
excess counts per Gaussian sigma of the smoothing function. 

No coiTection for the falloff of detection efficiency towards 
the edges of the field of view is applied because in the region 
of interest, around flie SNR RX JI713. 7-3946, the variation in 
gamma-ray acceptance is negUgible (see Fig. 1). 

2.5.2. Spectral analysis 

For the spectral analysis, the H.E.S.S. standard cut on the im- 
age size of 80 photo-electrons is applied to the data. To obtain 
a spectrum of a certain region in the sky, all events with re- 
constructed direction in that particular region are considered 
as ON events. A complication arises for the background es- 
timation. The gamma-ray acceptance and therefore the back- 
ground level depend strongly on energy; one cannot, as for the 
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image generation, simply use a 1-D radial lookup (which is in- 
tegrated over all energies) as a background estimate. Instead, 
the acceptance lookups would have to be generated in energy 
bins which in practice is difficult to handle. Another approach 
was applied here for the spectral analysis: background (OFF) 
events were selected from the same field of view, from the same 
data run, by selecting regions of the same size and form as the 
ON region, but displaced on a circle around the pointing di- 
rection of the system. The circle is chosen such that the OFF 
regions are at exactly the same offset (that is, at the same an- 
gular distance to the pointing direction) as the ON region. A 
minimum distance between the ON and OFF regions of 0.1° 
is required to avoid gamma-ray contamination. Furthermore, 
known gamma-ray sources in the field of view not associated 
with the test region are excluded from the OFF regions. This 
approach ensures that background events are taken at the same 
zenith and offset angles, which is crucial because of the de- 
pendence of the effective detection areas upon these two quan- 
tities, and it uses more or less the same region of sky, with 
similar night- sky-background noise. For an object of the size 
of RX J1713.7-3946 this results in one ON and OFF region, 
the latter being simply the reflection of the former at the system 
centre. This is illustrated in Fig. 5, where the OFF regions used 
for each observation position are drawn. 

After the geometrical selection, in order to obtain a differ- 
ential gamma-ray flux dN/dE in units of (TeV"^ cm~^ s~^), ON 
and OFF events (A^on and A^qff) are binned logarithmically in 
energy and divided by the mean effective area and the expo- 
sure time in each bin. The energy dependent effective area is 
determined for each data run, for the corresponding zenith (O^) 
and offset (6) angle range, multiphed by the live time of the 
run and added up. Then, for each energy bin i, the bin entries 
are divided by the width (Ai) of that bin to obtain a differen- 
tial flux value. The differential flux results from subtracting the 
differential OFF- from the ON-flux histogram: 

\ dE ji Ai 2runsrAi(ez, e) Ai 2runsrAi(0z, 9) 

The normalisation factor ai is determined from the ratio of the 
areas of the ON and OFF regions used during analysis of the 
whole data set and takes the exposure of the different observa- 
tions into account. 

An alternative approach is the event-by-event usage of the 
effective area, which served as a systematic check. Rather than 
determining a mean effective area for the whole data set, each 
event is weighted with the inverse of the effective area, taking 
the event zenith angle, offset and energy. 

A complication arises from the dependence of the effective 
detection area on zenith angle and offset. The Monte Carlo ef- 
fective areas are generated at certain discrete zenith angles and 
offsets. When generating the mean effective area corresponding 
to a certain zenith angle and offset range, the effective area is 
interpolated linearly between the simulation values (see Fig. 3, 
where the markers indicate values available from simulations 
and the lines are the linear interpolations). 
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Fig. 6. Correlation between the standard calibration, event reconstruc- 
tion, and background reduction method (Method 1) and an alternative, 
completely independent approach {Method 2, see main text for more 
details). The black solid and white dashed lines run through the ori- 
gin with a slope of 1 . A linear correlation is clearly visible, deviations 
from the line originate from different cut efficiencies of the two meth- 
ods. 



3. Results 

With the H.E.S.S. data set, the morphology of 
RX J 17 13. 7- 3 946 and its spectrum are resolved with 
high precision. Given that this is the first in depth analysis of 
such an extended source in VHE gamma rays, we present in 
the following first selected examples of extensive systematic 
tests that were performed in order to assure the stability of the 
analysis and then discuss the results. 

3.1. Morphology 

When analysing the morphology of an extended source, one 
aims for the best possible resolution with, at the same time, 
sufficient event statistics. In order to explore image structures 
and their stability, the data set was analysed using the same 
calibration and analysis software, but applying different sets 
of cuts (like accepting only three- and four-telescope events 
and events with image amphtudes larger than 80, 200, and 400 
photo-electrons), which resulted in different resolutions and 
different event statistics. Another important issue is the appro- 
priate modelUng and subtraction of the background. As a sys- 
tematic test, alternative background models have been applied 
and the results compared. First, a set of OFF runs, taken at the 
same zenith angles with very similar night- sky noise as the ON 
runs, has been used. Furthermore, the morphology was cross- 
checked using a completely independent caUbration and analy- 
sis approach, not only for the reconstruction of the shower ge- 
ometry but also for background estimation. Rather than using 
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Fig. 7. Gamma-ray image of RX J1713. 7-3946. The linear colour 
scale is in units of excess counts (see Sect. 2.5.1 for a description of 
image generation). The white contour lines indicate the significance of 
the different features, the levels are linearly spaced and correspond to 
5, 10, and IScr, respectively. The significance of each point has been 
calculated assuming a point source at that position, integrating events 
within a circle of 0.1° radius. In the lower left hand comer a simu- 
lated point source is shown as it would appear in this particular data 
set (taking the point-spread function and the smoothing into account) 
along with a black circle of 2' radius denoting the tr of the Gaussian 
the image is smoothed with. 



standard (Hillas) parameters for image parametrisation and re- 
constructing shower geometry based on these parameters, this 
approach is based on a 3-D modelling of Cherenkov photon 
emission during the shower development in the atmosphere as- 
suming rotational symmetry, thereby predicting pixel ampli- 
tudes (for more details, see Lemoine-Goumard & de Naurois 
2005). A comparison of the two analysis methods is shown in 
Fig. 6; plotted is the correlation of gamma-ray excess counts 
for the sky region around RX J1713. 7-3946. A linear correla- 
tion is clearly apparent illustrating the good agreement of the 
two independent methods. 

With the systematic tests mentioned above it could be 
shown that the main features of the gamma-ray morphology are 
stable when analysed with different cuts, different background 
models as well as with independently determined calibration 
coefficients and alternative analysis methods. 

Figure 7 shows a 2° x 2° field of view around 
RX J1713. 7-3946. A cut on the image size at 200 photo- 
electrons was applied resulting in a superior resolution of 0.08° 
(see Fig. 2). The corresponding energy range is w 300 GeV 
to 40 TeV. This image of RX J1713.7-3946 confirms with 
much higher statistics the 2003 H.E.S.S. measurement, shown 
for example in Fig. 1 of Aharonian et al. (2004b). There is no 
evidence for time variability, as expected for an object of the 
size of RX J1713. 7-3946. The overall gamma-ray appearance 




Fig. 8. Shown is a simple geometrical model for the emission from 
a thick sphere matched to the dimensions and relative fluxes of 
RX J1713. 7-3946. Left: 2-D projection of a thick and spherical 
radiating shell, 1° in diameter, smoothed with the H.E.S.S. point- 
spread function. Adapted empirically to match the radial shape of the 
H.E.S.S. data set, the dimensions of the geometrical sphere are 5.5 pc 
for the inner and 10 pc for the outer boundary if one assumes a distance 
of 1 kpc to the source, the presumed distance to RX J1713. 7-3946. 
The emissivity in the northern, western, and southwestern part is a 
factor of two higher than in the southeast and east. Right: Radial pro- 
file from the geometrical model compared to the H.E.S.S. data profile 
of RX J17 13.7-3946. The centre coordinates used for the data plot are 
ffj20oo = 17hl3m33.6s, (5j2ooo = -39°45'36". The geometrical model 
profile has been scaled to the same area as the data profile. 



resembles a shell morphology with bright emission regions in 
the western and northwestern part where the SNR is believed 
to impact molecular clouds (Fukui et al. 2003; Moriguchi et al. 
2005). It is worth noting that there is a possible gamma-ray 
void in the central-southeastern region. The cumulative signif- 
icance for the whole SNR is about 39cr with these hard cuts, 
which corresponds to an excess of * 7700 events from the re- 
gion of RX J1713. 7-3946. Drawn as white lines in Fig. 7 are 
in addition the contours of significance of the gamma-ray sig- 
nal (levels correspond to 5, 10, and 15cr). The significance has 
been calculated considering events that fall within an angle of 
0.1° of each trial source position. Thus, the contours quantify 
the significance for each point as if there was a point source at 
that position. The background estimate was derived from OFF 
runs as described in Sect. 2.5.1. The brightest parts of the SNR 
exceed 20cr. Except for the void structure, where the signifi- 
cance just exceeds 5cr, most of the remaining emission regions 
are well above lOcr. 

From the gamma-ray image presented here it is clear that 
the emission regions cannot be distributed homogeneously in 
the sphere RX J1713. 7-3946. The image is neither rotational 
symmetric nor does it exhibit a shallow peak towards the cen- 
tre. Instead, a shell seems to be apparent in the northern, and 
western to southwestern part. Apart from that, there is more or 
less uniform emission found in the rest of the SNR with a slight 
flux increase towards the southeastern boundary. This resem- 
bles very much the image one would expect from a thick spher- 
ical shell radiating gamma rays with enhanced emission on one 
side, as is illustrated in Fig. 8 where a geometrical model of a 
thick radiating sphere is presented. The good match in shape 
between the data and the toy-model profile lends support to 
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Fig. 9. Morphology of RX J17 13. 7-3946 as it appears at different energies. Shown from left to right are gamma-ray excess images with energies 
of £■ < 0.6 TeV, 0.6 TeV < E < 1.4 TeV, and 1.4 TeV < £. Drawn additionally as white lines are contours of significance, linearly spaced 
at 5, 10, 15cr (as in Fig. 7). Note the increase in the signal-to-noise ratio with increasing energy. The energy bands were chosen such that each 
band represents about a third of the full data set (taking events after cuts). Furthermore, all three images were smoothed with a Gaussian of 
2', which makes them directly comparable to each other, and to Fig. 7. The resolution in each energy band is indicated in the lower left hand 
comer of the images; the three data subsets have comparable resolutions of « 0.08° (the resolution of the intermediate energy band is about 
6% better). This might be counter-intuitive, given that at larger energies camera images get bigger and fluctuation efl"ects become negligible 
thereby improving the energy and direction resolution. However, in this case that effect is compensated by the increasing mean zenith angle of 
the large-energy events. 



the assumption that indeed it is the shell of RX J1713. 7-3946 
which radiates gamma rays. 

Figure 9 shows three images of RX J1713.7-3946 in three 
distinct energy bands, £■ < 0.6 TeV, 0.6 TeV < E < 1.4 TeV, 
and 1.4 TeV < E (left to right). The energy ranges were cho- 
sen such that each band represents a third of the data set. Note 
that the angular resolution of all three images is roughly the 
same which makes them readily comparable. The signal-to- 
noise ratio of the low-energy image is evidently smaller. The 
shell-like morphology of the SNR is slightly blurred by fluctu- 
ations. Correspondingly, the significance contours indicate that 
only the bright northwestern half is significant in this energy 
band. In contrast, the whole remnant sticks out significantly in 
the two higher-energy bands. Most of the northwestern parts 
exceed lOcr, the brightest spots even exceed 15cr for energies 
beyond 1 .4 TeV. 

From the visual impression the remnant does seem to emit 
gamma rays more uniformly with increasing energy. However, 
within errors, the radial shape appears to be the same in all three 
energy bands, as seen from the radial excess profiles around the 
centre of the SNR shown in Fig. 10. The morphology does not 
change significantly with energy. This is qualitatively compat- 
ible with the results of the spatially resolved spectral analysis 
(see section 3.2). 

3.2. Energy spectra 

The energy spectrum of RX J 17 13.7-3946 was measured with 
H.E.S.S. over a large energy range. Systematic tests on the 
shape and characteristics of the energy spectrum of the whole 
SNR included application of a slightly different spectral anal- 
ysis technique, different background models, analysis in dis- 
tinct data subsets like small and large zenith angles and the 
five observation positions, analysis applying different cuts on 
image intensity and telescope multiplicity, investigation of the 
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Fig. 10. Radial profiles around the centre of the SNR (q']2ooo = 
17hl3m33.6s, (5]2ooo = -39°45'36") in the three energy bands plotted 
in Fig. 9, generated from the raw (that is, not smoothed) gamma-ray 
excess images. The intermediate and high-energy band images have 
improved signal-to-noise ratios. The radial profiles from these energy 
bands have been scaled such that they have the same area as the low- 
energy profile (scaling factors are » 0.87 and a; 0.65, respectively), to 
enable direct comparisons. 



influence of the exact binning, of the energy estimation and the 
fit to the effective area histograms obtained from simulations. 
Furthermore, the results were cross-checked with the indepen- 
dent calibration and analysis scheme mentioned in Sect. 3.1. 
Representative examples for these systematic tests are shown 
in Fig. 11, where the spectrum of the whole SNR is plotted us- 
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Fig. 11. Shown are three spectra that were produced to explore the 
systematic uncertainties. The alternative spectra (blue squares and red 
triangles) were scaled by factors of 10"- and 10"'*, respectively, for 
presentation reasons. Upper set: Standard analysis and standard back- 
ground modelling using a background position from the same field of 
view, with the exact same shape and distance to the pointing direction 
as the signal region, but opposite to it, on the other side of the pointing 
direction. Middle set: Alternative spectral analysis technique (event- 
wise effective area weighting, see Sect. 2.5.2) with an independent 
background estimate taken from OFF runs. The background regions 
in the OFF data were again selected such that they have the same 
shape and distance to the pointing direction as the signal region in 
the RX J1713. 7-3946 observations. Lower set: Spectrum produced 
applying an independent analysis chain. The background was deter- 
mined similarly as for the upper set. A third spectral analysis tech- 
nique was applied here, described in Piron et al. (2001). Plotted as 
black line on top of all three spectra to guide the eye is the best fit of a 
power law with energy dependent photon index to the spectrum shown 
in the upper set (see Table 1 for details). The error bars on the spectral 
points denote ±lcr statistical errors. 



ing the standard background estimation from the same field of 
view, compared to a completely independent background es- 
timation based on OFF runs, and to the independent analysis 
chain (Lemoine-Goumard & de Naurois 2005). The three spec- 
tra are found to be well compatible with each other. 

In order to obtain a quantitative estimate for the system- 
atic error on each flux point, the background estimation, the 
spectral analysis technique (event-wise effective area weight- 
ing or average effective area determination per data run) and 
the absolute energy scale of the experiment were considered 
as dominant contributions. For the background modelling, an 
uncertainty of Aa = 1%, a being the background normalisa- 
tion factor, was derived from the RX J1713. 7-3946 data set 
by comparing the standard analysis to the analysis using an in- 
dependent background model derived from OFF runs, as men- 
tioned above. It should be noted that the 1% variation found 
in the total number of background counts is compatible with 
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Fig. 12. Differential gamma-ray energy spectrum of 
RX J1713.7-3946, for the whole region of the SNR (sohd black 
circles). The best fit of a power law with energy dependent photon 
index is plotted as black line. For comparison the H.E.S.S. 2003 data 
points are also shown (blue open circles). Note the vast increase 
in energy coverage due to the increased sensitivity of the complete 
telescope array. The spectrum ranges now from 190 GeV to 40 TeV, 
spanning more than two decades in energy. The data points reported 
by the CANGAROO-Il collaboration (Enomoto et al. 2002) for the 
northwestern part of the remnant are also shown as red triangles, the 
corresponding best fit result as dashed red line. Error bars are +ltr 
statistical errors. 

being due to a statistical variation and thus must be regarded 
as an upper limit. Taking the systematic uncertainty due to the 
background level and the two spectral analysis techniques, an 
energy dependent systematic eiTor was obtained by analysing 
the data set six times - scaling a by (1. + [-1,0, -t-1] x Ac) 
and applying both analysis techniques. The RMS of the result- 
ing six flux points in each energy bin was taken as the system- 
atic uncertainty. The uncertainty in the energy scale is a global 
uncertainty which might cause a shift of the whole spectrum 
to lower or larger energies. It is due to uncertainties in the at- 
mospheric transmission models used in simulations (see Funk 
et al. 2004) and uncertainties in the light collection efficiencies 
of the telescopes. The combined error is estimated to be 20%. 

When fitting a power law with index F to the spectrum, 
the systematic error on the integral flux obtained from the fit 
function is conservatively estimated to be 25%, on the fit index 
it is AF = 0.1. 

The spectrum of the whole SNR was determined by in- 
tegrating events within 0.65° radius around the centre of the 
SNR, anooo = 17hl3m33.6s, 5j2ooo = -39°45'36". OFF 
events were selected from a reflected region in the same field 
of view (see Sect. 2.5.2 for explanation and Fig. 5 for illustra- 
tion). To ensure optimum match in the offset distributions of 
ON and OFF events, runs taken directly on the source, where 
no appropriate OFF region can be selected in the same field of 
view, were discarded from the spectral analysis. Accordingly, 
the total live time reduced slightly to 30.5 hours. A size cut of 
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80 photo-electrons was applied for the spectral analysis. This 
results in a cumulative significance of 31cr coiTesponding to 
~ 15400 excess events (normalisation factor Q' = 1.3). 

The resulting spectrum of the whole SNR is shown in 
Fig. 12. The data is in excellent agreement with the previous 
measurement in 2003, which covered the energy range from 
1 TeV to 10 TeV. The latest data span more than two orders 
of magnitude in energy, from 190 GeV to 40 TeV. The best fit 
of a power law with energy dependent photon index is plotted 
(the exact formula is given below). It describes the data reason- 
ably well. Table 1 summarises fits of different spectral shapes 
to the data. Three alternative shapes have been used: a power 
law with an exponential cutoff E^, 
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In all cases, /q is the differential flux normalisation, the energies 
E are normalised at 1 TeV and photon indices are specified with 
F. 

All three alternative shapes describe the data significantly 
better than the pure power law. However, among the alternative 
spectral shapes, none is significantly favoured over the others. 
At the highest energies, above 10 TeV, there is still a significant 
gamma-ray flux in excess of 6cr. It should be noted, though, 
that in order to draw strong conclusions about the high-energy 
shape of the spectrum, more data with better statistics at the 
high-energy end of the spectrum are needed. 

The spectrum reported by the CANGAROO-II collabora- 
tion (Enomoto et al. 2002) for the northwest part of the SNR 
is also shown in Fig. 12. From a power-law fit to the data 
they quoted a photon index F = 2.84 + 0.15 (statistical) + 
0.20 (systematic) and a differential flux normalisation at 1 TeV 
/o = (1.63 +0.15 + 0.32) X 10-" cm-2 s"' TeV^^ The diff'erence 
between the two spectra is somewhat larger than the quoted er- 
rors of the measurements. However, the CANGAROO-II spec- 
trum is only for a part of the remnant. Moreover we note that 
the CANGAROO-II collaboration has recently revised their 
systematic eiTors upwards. For example, the Galactic Centre 
photon index, which was initially given as 4.6 + 0.5 (Tsuchiya 
et al. 2004), was recently quoted as 4.6;^j 2 (Katagiri et al. 
2005). 

Figure 13 illustrates the three spectral shapes that were 
found to describe the data reasonably well. The three curves 
are extrapolated to 1 GeV to compare them with the EGRET 
upper limit on the energy flux of 4.9 x 10"" erg cm"^ s 
ranging from 1 GeV to 10 GeV, centred at 2.9 GeV. The limit 
was determined at the H.E.S.S. position of RX J1713. 7-3946 
by modelling and subtracting the known EGRET source 
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Fig. 13. H.E.S.S. energy spectrum of RX J1713.7-3946. Plotted are 
the H.E.S.S. points with their ±1<t statistical errors in an energy flux 
diagram. The three curves (specified in Table 1) are the best fit re- 
sults of a power law with an exponential cutoff, a power law with 
energy dependent photon index, and a broken power law, extrapolated 
to 1 GeV to enable comparisons with the EGRET upper limit in the 
range of 1-10 GeV. The shaded grey band represents the systematic 
uncertainty on the measurement, originating from the uncertainty on 
the background estimation. The blue arrow indicates the 20% system- 
atic uncertainty on the energy scale, which might shift the whole curve 
in the given direction. 



3EG 1714-3857 (Hartman et al. 1999), assuming that 
3EG 1714-3857 is not linked to the gamma-ray emission of 
RX J1713.7-3946. Since flie H.E.S.S. location is in close 
vicinity of 3EG 1714-3857 (actually it is overlapping), this 
procedure could only be caiTied out successfully above 1 GeV. 
The systematic error band for the H.E.S.S. data was obtained as 
described above. It is centred on the mean value of the three fit 
curves and represents the systematic error due to background 
uncertainties only. The energy scale is an energy independent 
uncertainty; its scale and the direction the curve is shifted to 
are marked with a blue arrow at one representative position (at 
15 TeV). It is worth noting that the systematic uncertainty on 
the background has a considerable impact on the first few flux 
points because of the smaller signal-to-noise ratio (as compared 
to points at higher energies). For the spectrum shown here the 
systematic uncertainty is a; 18% for the two lowest-energy 
points; it decreases rapidly with increasing energy being well 
below 10% at 350 GeV. 

The results of the spatially resolved spectral analysis are 
shown in Fig. 14. It shows the distribution of photon indices 
over the SNR resulting from a power-law fit. The spectra 
were determined in rectangular regions, denoted 1-14, each 
0.26° X 0.26° in dimension. The fit range was limited to 8 TeV 
to account for (and avoid when fitting) the deviation from a 
power law seen in the spectrum of the whole remnant. Table 2 
summarises the fit results. There is a significant flux variation 
over the SNR. From the brightest region in the northwest to a 
relatively dim one in the central part, the flux varies by more 
than a factor of two. There is no significant difference in spec- 
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Table 1. Fit results for different spectral models. The differential flux normalisation Iq and the integral flux above 1 TeV (/(> 1 TeV)) are given 
in units of 10" cm -2 S-' TeV"' and 10"'^ cm ^ s respectively. The power-law fit is clearly an inappropriate description of the data, a power 
law with an exponential cutoff (row 2), a power law with an energy dependent photon index (row 3), and a broken power law (row 4; in the 
formula, the parameter S = 0.4 describes the sharpness of the transition from Ti to Ft and it is fixed in the fit) are equally likely descriptions of 
the H.E.S.S. data. Note that when fitting a broken power law to the data, some of the fit parameters are highly correlated. 

Fit Formula Fit Parameters (d.o.f.) /(> 1 TeV) 

IqE-^ /o = 17.1 ±0.5 F = 2.26 + 0.02 85.6 (23) 13.5 ±0.4 

loE-'^ exp(-E/E^) /o = 20.4 ±0.8 F= 1.98 ±0.05 £, = 12 ±2 27.4 (22) 15.5 ± 1.1 

£.-r+/j iog£ /o = 19.7 ±0.6 F = 2.08 ±0.04 yS = -0.30 ± 0.04 25.5 (22) 15.6 ± 0.7 

/ofi/f/' (i+£/£b'^) Io=OA+°-^2 F,= 2.06 ±0.05 F2 = 3.3 ± 0.5 £b = 6.7 ± 2.7 26.2(21) 15.4 ± 0.8 




Fig. 14. The image illustrates the results of the spatially resolved spectral analysis. Left part: Shown in red are gamma-ray excess contours 
from Fig. 7, linearly spaced at 30, 60, and 90 counts. Superimposed are the 14 boxes (each 0.26° x 0.26° in dimension) for which spectra were 
obtained independently. The dashed line is the 0.65° radius circle that was used to integrate events to produce a spectrum of the whole SNR. 
The photon index obtained from a power-law fit in each region is colour coded in bins of 0.1. The ranges of the fits to the spectra have been 
restricted to a maximum of 8 TeV (see Table 2 ). Right part: Plotted is the integral flux above 1 TeV against the photon index, for the 14 
regions the SNR was sub-divided in. The error bars are ±ltr statistical errors. Note that systematic errors of 25% on the flux and 0.1 on the 
photon index are to be assigned to each data point additionally. 



tral shape apparent, the photon indices agree with each other 
within statistical and systematic errors. The distribution of pho- 
ton indices has a mean value of 2.09 with a root-mean-square 
of 0.07. This is well compatible with the spectrum of the whole 
SNR when the fit range is also restricted to maximum 8 TeV 
for consistency (first row in Table 2). If one adds up the in- 
tegral fluxes above 1 TeV of the individual regions, a flux of 
15.1 X 10"'^ cm"^ s"' is obtained, 5% less than the flux of the 
whole SNR (with the restriction of the fit range). This is in 
excellent agreement with expectations; the boxes as they are 
plotted in Fig. 14 cover the region of RX J1713.7-3946 with 
significant gamma-ray excess almost completely. 

As can be seen from Fig. 14, right part, there is no correla- 
tion of the gamma-ray flux and the photon index visible in the 
data. This, together with the absence of any change in the spec- 
tral shape, is a remarkable difference between the gamma-ray 



and X-ray data. The spectral variation in X-rays was found to 
be much larger (see Cassam-Chenai et al. 2004). 

4. RX J1 71 3.7-3946 at other wavelengths 

4.1. Comparison with X-rays 

There is a striking similarity between the X-ray and the 
gamma-ray image of RX J1713. 7-3946, as they are shown in 
Fig. 15 for the 1-5 keV X-ray band. The overall morphology 
appears to be very similar, the brightest spots in both images 
are distributed on the shell, especially in the west. For a de- 
tailed comparison one must take into account the shghtly bet- 
ter resolution of ASCA compared to H.E.S.S.. For that pur- 
pose, the ASCA image was smoothed beyond the point-spread 
function of the instrument in order to match the H.E.S.S. res- 
olution. The smoothing radius was determined empirically by 
smoothing the ASCA point-spread function and comparing it 



F. Aharonian et al.: The y-ray supernova remnant RX J1713. 7-3946 



13 



• HESS 

-B-ASCA1-5keV 1 
-A- ASCA 5-10keV 



0.2 0.4 0.6 

_I 1 1 1 I 1 1 1 I 1 I 1 L 




0.2 



-a- > 



0.4 0.6 

J 1 1 1 1 L 



0.2 

_i 1 1 1 I I 



0.4 0.6 

1 I 1 1 1 L 




0.2 0.4 0.6 

_l 1 1 1 I 1 1 1 I 1 1 1 L 



-Q- 

^ 7 



0.2 0.4 0.6 

_l 1 1 1 I 1 1 1 I 1 1 1 L 



0.2 0.4 0.6 



8 



0.2 0.4 0.6 

_l 1 1 1 I 1 1 1 I 1 1 1 L 



r (deg) 

Fig. 16. Radial profiles for the eight regions marked in Fig. 15. Plotted are H.E.S.S. excess counts per unit solid angle (solid circles) as a 
function of distance r in degrees to the centre of the SNR, compared to soft (1-5 keV) and hard (5-10 keV) X-ray data. All data were corrected 
for relative acceptance. There is a very good general agreement between the keV and TeV data sets. The most pronounced differences appear 
in regions 4, where the TeV flux drops almost to zero at « 0. 15° from the centre, and 7, where a pronounced peak appears in the 1-5 keV X-ray 
data, which repeats neither in the 5-10 keV X-ray nor the TeV data. Note that the X-ray emission includes diffuse Galactic X-ray emission. 





Fig. 15. ASCA X-ray (1-5 keV band, Uchiyama 2005) image 
of RX J1713. 7-3946, overlaid with contours of the smoothed, 
acceptance-corrected H.E.S.S. gamma-ray image. The coloured con- 
tour levels are labelled and linearly spaced at 30, 60, and 90 counts. 
Drawn as gray thin lines are eight wedge-shaped regions for which the 
radial profiles are compared to each other in Fig. 16. Note that in the 
ASCA image, most of the regions (faint solid lines) do not reach as far 
as in the H.E.S.S. image (faint dashed lines), accounting for the lim- 
ited field of view of ASCA, whose coverage did not always extend to 
the boundaries of the SNR. As explained in the main text, the ASCA 
image was smoothed to match the H.E.S.S. point-spread function to 
enable comparison of the two images. 



Table 2. Fit results for distinct regions of the SNR. Given are for each 
region the photon index resulting from a power-law fit, the best-fit 
and the number of degrees of freedom (d.o.f.), the integral flux above 
1 TeV and the significance of the excess events in units of standard 
deviation (cr). The background for each region was determined from 
the same field of view, as described in Sect. 2.5.2, for each region 
separately, and hence the background estimates for different regions 
are not independent. The first row is the fit result of the whole SNR 
for comparison. For the whole table, the upper fit range was restricted 
to 8 TeV to avoid biases due to the deviation from a power law at high 
energies. 



Region 




r 


(d.o.f.) 


/(> 1 TeV) 


Excess 










(10-'-cm-2 s-') 


(0-) 


All 


2.12 


±0.03 


24.5 (18) 


15.9 ±0.6 


30.8 


1 


2.12 


±0.07 


34.7 (18) 


1.05 ±0.13 


12.9 


2 


2.24 


±0.06 


26.0(17) 


1.34 ±0.10 


17.2 


3 


2.11 


±0.06 


30.2(18) 


1.45 ±0.13 


16.7 


4 


2.10 


±0.08 


15.7(18) 


1.06 ±0.12 


11.5 


5 


2.09 


±0.07 


12.6(18) 


1.22 ±0.11 


13.3 


6 


2.13 


±0.06 


35.7(17) 


1.23 ±0.12 


14.1 


7 


1.95 


±0.08 


9.4(16) 


1.19 ±0.12 


10.9 


8 


2.11 


±0.12 


13.8(14) 


0.78 ±0.11 


8.0 


9 


2.11 


±0.10 


12.5 (16) 


0.89 ±0.11 


8.7 


10 


2.19 


±0.07 


24.8 (17) 


1.09 ±0.10 


14.1 


11 


2.08 


±0.08 


11.6(15) 


1.13 ±0.11 


11.8 


12 


2.01 


±0.11 


8.4(16) 


0.81 ±0.12 


8.2 


13 


1.98 


±0.10 


10.7(15) 


1.00 ±0.14 


9.8 


14 


2.08 


±0.11 


9.9(15) 


0.84 ±0.12 


9.4 



with the H.E.S.S. point-spread function. An optimum match 
was obtained for a smoothing radius of 0.037°. Furthermore, 
both images in Fig. 15 are corrected for relative detector accep- 
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tance. The X-ray-bright central point source IWGA J 17 13.4- 
3949, which was argued to be the compact relic of the super- 
nova progenitor (see for example Cassam-Chenai et al. 2004), 
was removed from the X-ray image for the purpose of compar- 
ison with gamma rays. 

After degradation of ASCA's resolution, the data were com- 
pared to each other in eight wedge-shaped regions as they are 
drawn in Fig. 15. In each wedge, radial profiles, that is, the 
number of counts per unit solid angle as function of distance 
to the centre, were determined. To account for the differences 
in the absolute count level the X-ray images were scaled by a 
normalisation factor, which has been calculated as the ratio of 
TeV counts, integrated in a rectangle encompassing the SNR 
(and within the ASCA field of view), to keV counts, integrated 
in the same rectangle. The result is shown in Fig. 16 for all 
eight wedges. The overall good agreement in shape of the pro- 
files is clearly visible, differences appear only at a few places, 
for example in region 4 and 7. For a quantitative statement on 
the compatibility of the two data sets, however, one would have 
to model and subtract the contribution from Galactic diffuse X- 
ray emission in the ASCA image, which amounts presumably 
to 10% or less in the X-ray bright parts of the SNR, but might 
increase to 30% in the faint parts in the east. 

The interesting question of the boundaries of the SNR and 
if they are the same in X-rays and gamma rays can unfortu- 
nately not be addressed with the ASCA data set due to limited 
sky coverage. 

4.2. CO and radio observations 

CO data at 2.6 mm wavelength of RX J1713.7-3946 and its 
surroundings were taken with the 4-m, mm and sub-mm tele- 
scope NANTEN in 2003 (Fukui et al. 2003). Based on posi- 
tional coincidences of CO and X-ray peaks and (in velocity 
space) shifted CO peaks, Fukui et al. (2003) concluded that 
the SNR blast wave is interacting with molecular clouds sit- 
uated on its western side at a distance of 1 kpc. Further 
possible support for this scenario was recently published in 
Moriguchi et al. (2005), where high gas excitations are reported 
for this part, which could arise from heating of the molecular 
gas by the shock wave. The CO intensity distribution in the 
corresponding velocity interval is shown in Fig. 17, together 
with H.E.S.S. gamma-ray excess contours. One notes that in 
the central and southeastern part of the SNR the CO emission 
becomes very faint or is completely absent. Apart from that, 
there are local CO maxima that coincide with TeV-bright parts 
on the western side of the SNR. The azimuth profile plotted 
on the right-hand side of Fig. 17 illustrates a global agreement 
between the two measurements, regions with low gamma-ray 
flux reveal also low CO intensity, but there is no exact propor- 
tionality between the two measurements for the shell region of 
RX J1713. 7-3946. Taking the peak values, one notes that they 
are shifted with respect to each other and that the gamma-ray 
flux varies by a factor of about three, whereas the CO intensity 
drops by roughly a factor of 100 in the central-southeastern 
part. 




17h16m 17h14m 17h12ni 
Right Ascension (hours) 

Fig. 18. 1.4 GHz ATCA radio image (Lazendic et al. 2004, courtesy 
of P.Slane). The linear colour scale is in units of Jy beam"'. Overlaid 
are coloured contour lines of the H.E.S.S. gamma-ray excess image. 

Figure 18 shows a comparison of the 1.4 GHz radio image 
obtained with ATCA (Lazendic et al. 2004) and the H.E.S.S. 
gamma-ray excess contours. The SNR is very faint in the radio 
band, there are two faint arc-like structures of emission to the 
west of RX J 17 13.7-3946, almost perpendicular to each other, 
one of them directly coincident with the brightest TeV emission 
region. There is no notable resemblance between the two wave- 
length regimes. Spectral analysis of the X-ray- and TeV-bright 
northwestern part of the SNR shell yields a spectral index of 
0.50 + 0.40, derived from two flux values taken at 1.4 GHz and 
2.5 GHz. This measurement is used further down when com- 
paring the H.E.S.S. spectral data to broadband models. 

5. Possible emission processes 

One of the key issues in the interpretation of the observed 
gamma-ray emission is the identification of the particle pop- 
ulation responsible for the generation of the gamma rays. The 
close correlation between X-rays and gamma rays might in- 
dicate an electronic origin; models of supernova remnants as 
Galactic cosmic -ray sources, on the other hand, suggest that 
primarily a hadronic component from pion decays exists. To 
identify the different contributions, the wide-band electromag- 
netic spectra from radio to multi-TeV gamma-rays must be 
compared to model calculations. 

In the literature, different schemes are employed to model 
broadband emission from SNRs. Phenomenologically oriented 
models (Mastichiadis & de Jager 1996; Aharonian & Atoyan 
1999) start by ad hoc assuming particle acceleration spectra - 
usually as power laws with a cutoff - to derive particle spectra 
taking into account energy losses and then calculate the elec- 
tromagnetic spectrum with additional assumptions concerning 
the local magnetic field, the radiation fields which serve as tar- 
get for the IC process, and the gas density. Spectral parame- 
ters are either taken from acceleration models, resulting in a 
spectral index around 2, or determined from data. More so- 
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Galactic Longitude (deg.) Azimuth (deg.) 

Fig. 17. Left panel: Shown are the intensity distribution of CO {J = 1-0) emission (Fukui et al. 2003) (linear colour scale in units of K km s"' , 
truncated at a value of 23 to highlight important features), derived by integrating the CO spectra in the velocity range from -11 km s"' to 
-3 km s"' (which corresponds to 0.4 kpc to 1.5 kpc in space). Overlaid are coloured contours of the H.E.S.S. gamma-ray excess image. The 
levels are labelled and linearly spaced at 30, 60, and 90 counts. Note that the image is shown in Galactic coordinates. Right panel: Azimuth 
profile plot, that is, number of counts as a function of the azimuthal angle, integrated in a 0.2°-wide ring covering the shell of RX J1713. 7-3946 
(dashed yellow circle in the left-hand panel). Plotted are the H.E.S.S. gamma-ray and the NANTEN data set. 



phisticated gamma-ray models account for the non-linear ef- 
fects arising from the interaction of the accelerated particles 
with the shocked supernova shell, which result in deviations 
from pure power laws, with spectra flattening at higher ener- 
gies (Berezhko & Volk 1997; Baring et al. 1999). 

The original discovery paper of VHE gamma-ray emis- 
sion from RX J1713.7-3946 (Muraishi et al. 2000) claimed 
electrons as the likely source particle population. However, it 
soon became evident that a consistent modelling of the spec- 
tra is hard to achieve in simple one-zone models. Apart from 
the choice of the electron spectrum, the only free parameter is 
the magnetic field strength, which controls the spacing of the 
synchrotron and IC peaks in the SED together with their rela- 
tive intensities and - one should add - the amount of radiative 
cooling of the accelerated component. Enomoto et al. (2002) 
noted that for modest magnetic fields - B equal to a few yuG 
- the measured intensity ratios are reproduced but the gamma- 
ray spectra are much too hard. Using higher fields, one can 
match the gamma-ray spectra at the expense of dramatically 
increased X-ray yields. While the H.E.S.S. data differ from the 
CANGAROO-II data both in terms of the region covered and 
the exact values for flux and index, this conclusion for the elec- 
tronic scenario remains basically valid. The agreement can be 
improved by introducing an additional parameter to decouple 
the X-ray intensity and the spectral shape, namely the mag- 
netic field filling factor which allows the X-ray flux to be tuned 
without change of the spectra. With very small filling factors of 
0.001 (Pannuti et al. 2003) to 0.01 (Lazendic et al. 2004), dif- 
ficult to justify physically, the X-ray and the CANGAROO-11 
gamma-ray spectra can be described for magnetic fields around 
10 fiG to 15 juG in the emitting region. This latter approach is 
not followed here. 



The validity of electronic models could be judged more eas- 
ily if the magnetic field values in the remnant were known. For 
typical shock compression ratios around 4 and pre-shock inter- 
stellar fields of a few fiG, fields of 10 fiG to 15 fiG are at the 
lower limit of the expected range; mechanisms of dynamical 
field amplification in non-linear shocks (Lucek & Bell 2000; 
Befl & Lucek 2001; Bell 2004) will generally result in higher 
fields. The narrow filaments visible in many high-resolution X- 
ray images of SNRs (see, e.g., Bamba et al. 2005) have been 
pointed out to provide means to probe magnetic fields (Vink 
& Laming 2003; Berezhko et al. 2003): only relatively high 
fields can result in sufficiently rapid cooling of electrons to 
make such filamentary features possible and visible. On the ba- 
sis of the structures seen in Chandra images in the northwest of 
RX J1713.7-3946 (Uchiyama et al. 2003), Volk et al. (2005) 
have argued that fields between 58 fiG and a few 100 fiG might 
be possible, depending on the detailed assumptions about the 
remnant's morphology'. Such high fields - likely to be present 
throughout the remnant - would rule out a leptonic origin of 
VHE gamma rays right away. 

On the basis of the difficulty of accommodating broad- 
band spectra in a single-zone electronic model, Enomoto 
et al. (2002) proposed RX J1713.7-3946 as the first well- 
identified proton accelerator. This interpretation was criti- 
cised by Butt et al. (2002) and Reimer & Pohl (2002) since 
the CANGAROO-II spectra, extrapolated to lower energies, 
would violate the flux level of the nearby EGRET source 
3EG 1714-3857 (Hartman et al. 1999), which, if not associ- 

' We note that in Hiraga et al. (2005), Fig. 2, a radial profile from 
the XMM image is shown. It reveals another very thin filament-like 
structure in the west of RX J1713. 7-3946 which is a sign of high 
magnetic field values. 
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Fig. 19. Broadband SED of RX J1713.7-3946. The ATCA radio data 
and ASCA X-ray data (ffiraga 2005) for the whole SNR are indicated, 
along with the H.E.S.S. measurement and the EGRET upper limit. 
Note that the radio flux was determined in Lazendic et al. (2004) for 
the northwest part of the shell only and was scaled up by a factor of 
two here to account for the whole SNR. The synchrotron and IC spec- 
tra were modelled assuming a source distance of 1 kpc, an age T of 
1000 years, a density « of 1 cm"^, and a production rate of relativistic 
electrons by the acceleration mechanism in the form of a power law 
of index a = 2 and an exponential cutoff of £0 = 100 TeV. Shown 
are three curves for three values of the mean magnetic field: 7 yuG, 
9 yuG, and 11 yuG, to demonstrate the required range of the B field 
strength for this scenario. The electron luminosity is adopted such that 
the observed X-ray flux level is well matched. For the three magnetic 
field values the luminosity Le is Le = 1.77 x 10'' erg s"' (7 yuG), 
Le = 1.14 X 10" ergs-' (9 yuG), and Le = 0.81 x 10" erg s"' (UyuG). 

ated with RX J1713. 7-3946, must then be considered as upper 
limit on the GeV emission. However, the EGRET limit can be 
circumvented by reducing the amount of low-energy protons 
compared to the E"^ spectrum. This can be achieved by the ad 
hoc assumption of a spectral break, or - for the CANGAROO- 
II data - by assuming a flatter overall spectrum with a photon 
index smaller than 2. 

A very detailed modelling is beyond the scope of this paper; 
the models presented in the following serve mainly to illustrate 
that spectra and energetics can be reproduced with plausible 
input parameters. 

5A. Electronic scenario 

In Fig. 19 the synchrotron and IC emission from relativistic 
electrons are modelled within the framework of a one-zone 
model in which the electron acceleration and gamma-ray emis- 
sion take place in the same region. It is assumed that the pri- 
mary electrons follow a power law with index a = 2 and with 
an exponential cutoff Eq, 

e(£) = eo£-"exp(-£/£o), 

and that they are produced continuously over a fixed time in- 
terval T inside a region with given homogeneous distributions 
of magnetic field strength B and ambient gas density n. The 
energy distribution of the electrons is then calculated taking 
into account energy losses due to IC and synchrotron emission, 
Bremsstrahlung and ionisation as well as losses due to Bohm 



diffusion. The broadband energy distribution of the source is 
calculated for an age ofT- 1000 years, an average gas density 
of n = 1 cm"^, and a distance to the source of D = 1 kpc. For 
the IC emission, canonical interstellar values for the seed pho- 
ton densities were considered: Wcmb = 0.25 eV cm"^ for the 
cosmic microwave background (CMB), Wsl = 0.5 eV cm"-' 
for optical star light and Wir = 0.05 eV cm"^ for infrared 
background light. The absolute electron production rate, Qo, 
is determined from the constraint of matching the observed X- 
ray flux level. Figure 19 shows the resulting model curves, to- 
gether with measurements in various wavelength regimes, for 
three different average magnetic field values. From the abso- 
lute levels it is evident that a magnetic field around 10 yuG is 
required in order to explain both the X-ray and gamma-ray flux 
levels. On the other hand one notes that such a model with the 
above mentioned parameters does not provide a reasonable de- 
scription of the H.E.S.S. data. The IC peak appears too narrow 
to reproduce the flat TeV emission. The detailed inclusion of 
non-linear acceleration effects should not change the situation 
very much. They are expected to steepen the synchrotron SED 
above the radio range. Synchrotron cooling of the accelerated 
electrons then tends to produce a flat-topped synchrotron and 
accordingly IC maximum. It is, however, a long way to flatten 
the IC spectrum so extensively at low energies as to achieve 
agreement with the H.E.S.S. spectrum. 

Obviously, the simple model presented here served basi- 
cally to underline the main arguments. Nevertheless, the con- 
clusion that a power-law production spectrum fails to simulta- 
neously account for the radio. X-ray and gamma-ray data ap- 
pears to be a generic and stable feature; additional parameters 
are required to decouple either the TeV and X-ray/radio fluxes 
- such as a fillingf actor - or the X-ray and radio spectra - such 
as an ad-hoc spectral break, which for the given source age and 
magnetic field can not be justified as an effect of radiative cool- 
ing. 

5.2. Hadronic scenario 

Assuming alternatively that nuclear cosmic-ray particles, ac- 
celerated at the SNR shock, dominantly produce VHE gamma 
rays, theoretically the most plausible differential energy spec- 
trum of accelerated nuclei is a concave £'"'^*^^-type spectrum, 
due to nonlinear back coupling, with a cutoff at gamma-ray en- 
ergy Ec, where F is decreasing towards higher energies (just 
below the TeV energy range) to a value between 1 .5 and 2, be- 
fore the spectrum is steepening again in the cutoff region. In the 
test-particle approximation one expects F ^ 2. The H.E.S.S. 
spectrum is indeed compatible with such a scenario. Figure 20 
shows a vFy representation of the H.E.S.S. data, together with 
the best-fit curve of a power law with an exponential cutoff (see 
Sect. 3.2, Table 1), extrapolated to small energies. Compared 
to that a curve is plotted which takes the kinematics of the 
production process of gamma rays, pp li' — » yy, into ac- 
count. The power law spectrum continues to smaller energies 
with an index of ^ 2, as expected in the test-particle limit, until 
the suppression of gamma rays due to Tr^-decay kinematics is 
encountered and the curve is turning down. Note that already 
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Fig. 20. H.E.S.S. data points plotted in an energy flux diagram. Tliey 
shaded grey band is the systematic error band for this measurement 
(see Sect. 3.2). The black curve is the best fit of a power law with ex- 
ponential cutoff to the data, extrapolated to lower energies. The dashed 
blue curves is the same function, but it takes the tt" kinematics into ac- 
count. The EGRET upper limit from 1 GeV to 10 GeV is plotted as 
red arrow. 

the extrapolation of the H.E.S.S. spectrum is well below the 
EGRET upper limit from the position of RX J1713.7-3946, 
introduced in Sect. 3.2. Taking into account non-linear effects 
would harden the gamma-ray spectrum even more. 

One should mention at this point that on the theory side 
other mechanisms to suppress contributions from low-energy 
(E around 10 GeV) protons have been considered (e.g., Malkov 
et al. 2005). When the particles upstream of the shock hit a 
dense target with a spatial gradient, such as a dense molecular 
cloud, the gamma-ray emission of low-energy protons might 
also be suppressed due to the energy dependence of the diffu- 
sion length. In a more general context, such mechanisms - an 
accelerator of finite lifetime interacting with a target at a dis- 
tance where diffusion time scales are comparable to the source 
lifetime - have been studied by Aharonian & Atoyan (1996). 
Such arguments, however, need to be reassessed when the exact 
location of the clouds with respect to the shock front is known. 

To calculate the energetics in a hadronic scenario, the mean 
target gas density available for gamma-ray production in the 
region of RX J1713. 7-3946 is a key question. The CO image 
shown in Fig. 17, left-hand side, reveals a hole in the CO emis- 
sion and accordingly in the molecular hydrogen distribution in 
the central and eastern part of the SNR. In contrast, the TeV 
emission fills the whole region of the SNR (see Fig. 7). As is 
shown on the right-hand side of Fig. 17, there is no exact corre- 
lation between VHE gamma rays and CO intensity. From this 
one can conclude that in all likelihood cosmic rays do not pen- 
etrate the clouds uniformly. The bulk of the VHE gamma rays 
is not linked to the molecular clouds but must be due to inter- 
actions with a different target. Indeed, the rather good spheri- 
cal shape of the remnant together with the fact that the X-ray 
and gamma-ray emission only varies by a factor of two to four 
across the remnant lends further support to the scenario that the 



SNR is running into a more or less uniform and probably low- 
density environment. Although it seems to be beginning to in- 
teract with the dense clouds to the west, the ones that are seen 
by NANTEN, VHE gamma rays are dominantly produced in 
cosmic-ray interactions with rather uniform ambient gas. One 
possible scenario is that the SNR is the result of a core-collapse 
supernova explosion that occurred into a wind bubble of a mas- 
sive progenitor star An SNR shock expanding into the bubble, 
with an ambient density much lower than that suggested by 
an average molecular cloud scenario, could explain the rela- 
tive uniformity of the gamma-ray emission, compared to the 
large density variations in the clouds which likely surround the 
remnant (for a theoretical treatment of such configurations, see 
Berezhko & Volk (2000)). 

The local target density is a crucial parameter in this sce- 
nario. With the NANTEN measurement of the void in the cen- 
tral part of the SNR one might constrain the local density in 
that region. The sensitivity of the final NANTEN data set as 
quoted in Fukui et al. (2003) corresponds to a molecular col- 
umn density of 8.3 x 10'^ cm"^ assuming the conventional 
conversion relation from CO intensity to H2 column density 
(0.4 K km s"'). Taking the diameter of the SNR as w 20 pc for 
1 kpc distance, one can deduce an upper limit on the molec- 
ular hydrogen column density of x 2.6 cm"^ in parts of the 
SNR without detectable CO emission. The other existing con- 
straint was inferred in Cassam-Chenai et al. (2004) from XMM 
data, based on the lack of thermal X-ray emission. By fitting the 
spectra with an absorbed power-law model and adding a ther- 
mal component, an upper limit on the mean gas temperature 
and, important here, the mean hydrogen number density of the 
ambient pre-shock medium of 0.02 cm"^ was obtained. One 
should note, though, that this value is likely to be too low - if 
the shocks are strongly modified by the accelerated particles, 
the shock heating is substantially reduced and the data would 
be consistent with higher densities. 

Assuming for now a mean target gas density of « ^ 1 cm"^, 
uniformly spread throughout the remnant, in accordance with 
the NANTEN, but not the XMM limit, one can calculate the 
proton energetics implied by the gamma-ray flux measured 
from 0.2 to 40 TeV. The total energy in accelerated protons 
from about 2^00 TeV, required to provide the observed flux, 
can be estimated as Wj,°'(2^00 TeV) * tj,p^„o x Ly(0.2- 
40 TeV), where tpp^^ w 4.5 x 10'^(«/1 cm"^)"' s is the char- 
acteristic cooling time of protons through the jfi production 
channel, Ly(0.2^0 TeV) = 4Kd^Wy(0.2-40 TeV) is flie lumi- 
nosity of the source in gamma rays between 0.2 and 40 TeV, 
and Wy{Q.2-4Q TeV) is the gamma-ray energy flux for the cor- 
responding energy range. Assuming then that the proton spec- 
trum with spectral index o- « F continues down to 1 GeV, the 
total energy in protons can be estimated and compared to the 
total assumed mechanical explosion energy of the supernova of 
10^' erg. These calculations reveal very similar values for the 
three spectral shapes given in Fig. 20 in the gamma-ray energy 
range between 0.2 and 40 TeV: for the gamma-ray energy flux 
one obtains Wy(0.2^0TeV) ^ 10"'° erg cm"^ s"', the gamma- 
ray luminosity is Ly(0.2-40 TeV) * 10^^ (^j^A—f erg s"', 
and the corresponding energy content of protons is Wp°'(2- 
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400 TeV) « 6 x lO^"^ {j^f (^) ' erg. The resulting 
total energy in protons, after extrapolating the proton spectrum 

to 1 GeV and using £51 = 10^' (to^)' (t^)"' erg, is 
then W]f' «i 0.19 X £51 for a power law with exponential cutoff, 
W^"' * 0.08 X £51 for a power law with energy dependent index, 
and W^°' « 0.26 x £51 for a broken power law. These numbers 
are consistent with the notion of an SNR origin of Galactic 
cosmic rays involving the canonical * 10% conversion effi- 
ciency of the total supernova explosion energy. The H.E.S.S. 
gamma-ray flux level is close to what was predicted in Drury 
et al. (1994) from nearby young SNRs for ambient densities of 
n « 1 cm"^. One should keep in mind though that the order of 
magnitude uncertainties in the measurements of the distance to 
the source d and of the local gas density n feed directly into 
these estimates. 

5.3. Discussion and conclusions 

The models and ideas presented in this section were aiming at 
exploring the possibilities available in explaining the observed 
VHE emission in purely electronic and purely hadronic scenar- 
ios. It is found that in the hadronic scenario, assuming gamma 
rays to stem from decays, the extrapolation of the H.E.S.S. 
spectrum to lower gamma-ray energies leads to a picture that 
is consistent with the low-energy EGRET data. Furthermore, 
the spectral shape is well compatible with cosmic-ray acceler- 
ation theory. The energy requirements implied by the gamma- 
ray flux are in agreement with expectations from cosmic-ray 
acceleration in shell-type SNRs in our Galaxy, if one assumes 
a local target gas density of n « 1 cm"^ and takes the cur- 
rently preferred distance estimate of 1 kpc. Unfortunately both 
of these parameters are not very well measured. The distance 
estimate, which factors quadratically into the energetics cal- 
culation, has uncertainties in the order of at least 30%. For 
the local target density there exist only upper limits, since 
from comparisons with CO data it turns out that gamma rays 
are most likely not exclusively linked to the dense molecular 
clouds surrounding the SNR. These clouds, however, obscure 
the measurement of the actual local target material available 
for gamma-ray production. Only towards the interior and the 
southeast of the SNR, where there is a hole in the molecular 
colunm density, is there hope to actually measure and con- 
strain the density. Existing estimates in these regions are the 
NANTEN upper hmit of 2.6 cm"^, which does not cause any 
problem with the assumption made above, and the XMM up- 
per limit of 0.02 cm"^ which, if correct, would seriously chal- 
lenge the idea of a hadronic scenario of gamma-ray production 
at least for this object. 

In the electronic scenario, on the other hand, the data are 
not easily reproduced taking only IC emission into account. 
The very low magnetic field oi ^ 10 juG, fixed by the ratio 
of synchrotron to IC flux, exceeds typical interstellar values 
only sUghtly and is difficult to reconcile with the paradigm of 
the difliisive shock acceleration of cosmic rays at supernova 
shock waves which predicts strong field amplifications in the 
region of the shock (Lucek & BeU 2000; BeU & Lucek 2001; 
Bell 2004). In the case of RX J1713.7-3946 it was indeed con- 



sidered possible by Volk et al. (2005) that the magnetic field 
strength at the SNR shock front significantly exceeds typical 
interstellar values. 

Complete understanding of gamma-ray emission processes 
can only be achieved by taking a broadband approach and us- 
ing all the available measurements in the difi^erent wavelength 
regimes. In Sect. 4 the TeV data set was compared to X-ray, 
radio and CO emission measurements of the region surround- 
ing RX J1713.7-3946. While there is no obvious resemblance 
with the radio image, it turns out that there is a striking spatial 
correlation between the ASCA X-ray and the H.E.S.S. gamma- 
ray data. Most of ffie emission regions seem to exhibit exactly 
the same morphology in boffi wavelength regimes. At first sight 
this supports the idea that X-rays and gamma rays are produced 
by the same particle population, namely electrons. Assuming a 
constant magnetic field throughout ffie remnant (not the most 
likely configuration), the intensity (and spectrum) of both syn- 
chrotron and IC radiation trace the density (and the spectrum) 
of electrons, giving rise to ffie observed correlation. If the VHE 
gamma rays were due to non-ffiermal Bremsstrahlung of elec- 
trons, which is correlated with gas density, ffie observed corre- 
lation could be due to a magnetic field and gas density correla- 
tion. However, as can be seen from Fig. 19, Bremsstrahlung 
dominates over IC radiation only for very large values of 
nn > 100 cm"^, which are not compatible wiffi the CO mea- 
surements from ffie centre of the SNR, as mentioned above. But 
even given such a high density it is questionable if density, field 
strength and electron spectra can be fine-tuned such that the ex- 
perimental results are approximately reproduced. Anoffier dif- 
ficulty for an electronic interpretation arises from ffie observa- 
tion by Cassam-Chenai et al. (2004) that ffie X-ray spectra are 
steeper in ffie presumed shock front in the west, where ffie blast 
wave probably impacts ffie molecular cloud, than in the souffi- 
east, where ffie front propagates into a low density medium. It 
is very difficult to explain why the spectral shape in X-rays, but 
not in gamma rays, changes significantly in distinct regions of 
the shock, if ffiey stem from the same particle population. If 
on ffie offier hand gamma rays originate dominantly from nu- 
cleonic cosmic rays, a spatial correlation between X-rays and 
gamma rays is not automatically ensured eiffier. There are two 
possible scenarios. The correlation could point to a common 
acceleration process accelerating both electrons and protons - 
indeed expected in the theory of diffusive shock acceleration - 
such ffiat ffie spatial distributions are to first order ffie same and 
only differ because of the different loss processes. The second 
alternative is a correlated enhancement of magnetic field and 
local gas density. 

Anoffier possibiUty of course is that ffie VHE gamma rays 
are a roughly equal mixture of two components, produced by 
both electrons and protons. However, this scenario seems un- 
likely since the energy-independent gamma-ray morphology 
and ffie absence of variations in spectral shape would again re- 
quire fine-tuning of parameters like ffie magnetic field B and 
the ambient density n. 

We conclude that ffie straightforward and simplest ap- 
proaches in both scenarios lead to problems and one has diffi- 
culties in finding unequivocal evidence for eiffier of them when 
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using all the available broadband data. Nevertheless, the shape 
of the gamma-ray spectrum favours a hadronic scenario. 

6. Summary 

The VHE gamma-ray emission of RX J1713.7-3946 was mea- 
sured with unprecedented precision with H.E.S.S.. The accu- 
racy of the measurement is now approaching the level of X-ray 
measurements of this source, with the distinct advantage that 
H.E.S.S. covers the whole SNR within its field of view. With 
the 2004 data, a close spatial correlation between X-rays and 
gaimna rays was found. The overall gamma-ray energy spec- 
trum was measured over more than two decades. There are in- 
dications for a deviation from a pure power-law spectrum. The 
data seem to be reasonably well described by a power law with 
an exponential cutoff and a power law with energy dependent 
photon index, as well as a broken power law. At the current 
stage further investigations about the shape of the spectrum at 
the highest energies accessible to H.E.S.S. are hampered by 
the limited event statistics. The large data set has allowed for 
a spatially resolved spectral study. No significant variation in 
the gamma-ray spectral shape over the SNR is found. The flux 
varies by more than a factor of two across the SNR. The north- 
em and western parts of the shell, where the SNR is believed 
to impact molecular clouds, are significantly brighter than the 
remaining parts. 

We presented broadband modelUng ideas and discussed 
RX J1713.7-3946 in terms of the available data from all wave- 
length bands including the H.E.S.S. gamma-ray signal. Two 
scenarios were addressed, one where gamma rays originate 
from electrons and one where they originate from protons. In 
both cases the large uncertainties on crucial parameters like the 
magnetic field strength and the effective ambient density, which 
are not directly accessible to measurements, hamper decisive 
conclusions. Nevertheless, the proton scenario is favoured be- 
cause of the shape of the gamma-ray spectrum. From the theory 
side, the remaining challenge is the connection of the differ- 
ent particle species, VHE electrons and nuclei, in a consistent 
broadband model of RX J1713.7-3946. Experimentally, with 
the current gamma-ray data set, more precise measurements of 
the surrounding molecular clouds are clearly needed in order to 
link emission regions of VHE gamma rays to regions of known 
density. 
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